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In the GRID technique one measures 
Doppler-broadened line profiles of 7 transi- 
tions using the high resolution crystal 
spectrometers GAMS, which are installed 
at the high flux reactor of the ILL Greno- 
ble. One of the essential applications of 
this technique is the measurement of 
nuclear state lifetimes. In the present 
contribution the precision and the principal 
limits of the technique are discussed. 
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1. Introduction 

During the last three decades, crystal spectrometers 
have been developed and installed at the high flux reac- 
tor of the Institut Laue-Langevin. Currently, three 
devices are available: the single bent crystal spectrome- 
ter GAMS 2/3 and the two double flat crystal spectrom- 
eters GAMS4 and GAMS5. They are installed in line to 
an in-pile tangential beam tube. In this beam tube 
targets can be placed at a distance of 55 cm to the 
reactor core where the thermal neutron flux equals 
5 X lO''' s"' cm"^. A general scheme of the experimental 
set-up is shown in Fig. 1. The combination of the 
intense 7 flux and the ultra high resolution of the double 
flat crystal spectrometers GAMS4 and GAMS5 has al- 
lowed the development of the Gamma Ray Induced 
Doppler broadening technique. Currently, this 
technique is widely used for nuclear state lifetime 
measurements and the study of atomic collision pro- 
cesses in solids [2]. It is the purpose of this contribution 
to give a brief description of the technique itself. 
Furthermore, an attempt to estimate quantitively the 
sensitivity of the technique in nuclear state lifetime 
measurements is presented. 



2. GRID Basics 

A detailed description of the GRID technique has 
been given in a series of publications [1,2,3] and there- 
fore only the basic ideas should be presented here. 

Thermal neutron capture produces an excited 
nucleus, which deexcites either by 7 emission or other 
decay processes. Whenever an emission process takes 
place, the nucleus is recoiling due to momentum con- 
servation. In the particular case of a 7 cascade this 
means that the measured energy of any secondary 7 ray 
might be Doppler shifted by a recoil induced during a 
prior emission of another 7 ray. Due to the spatial distri- 
bution of the recoil directions, it is a Doppler broaden- 
ing, which will be obtained. The broadening depends on 
the velocities of the recoiling nuclei at the time of the 
secondary emission and therefore on three factors: i) the 
initial recoil velocity, ii) the time between the initial 
recoil and the emission of the observed 7 ray e.g., the 
lifetime of the intermediate level and iii) the slowing 
down of the recoiling particles. The measurement of the 
Doppler-broadening allows extraction of information 
on one of these quantities, assuming the other two to be 
known. 
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Fig. 1. Experimental set-up of the crystal spectrometers installed at the in-pile tangential beam tube of the high flux reactor of the ILL. 



For illustration, we consider first the very simple case 
of a two step y cascade, where the intermediate level is 
fed to 100 % by a primary transition of energy £.,1. As 
the initial recoil velocity is known to be Vr = Eyi/Mc, 
where M is the mass of the recoiling isotope and c the 
velocity of light, the Doppler-broadened lineshape of 
the secondary transition Ey2 depends only on the 
lifetime of the intermediate nuclear state and the 
slowing down of the recoiling atoms. 

In order to have a description for the slowing down 
process of the recoiling atoms Jolie et al. [2, 4] have 
assumed a model, in which the slowing down process is 
assumed to be isotropic. Here, the description of the 
recoil motion is reduced to the determination of the time 
dependence of the average velocity v{t) of all recoiling 
atoms. This assumption is strongly based on the use of 
powder targets. In this case, the y rays are emitted from 
a large number of microcrystals and any microscopic 
anisotropy of the slowing down process is averaged out 
due to the random alignment of the microcrystals. 
Assuming v{t) to be known, the Doppler-broadened y 
lineshape is given by the equation 



Io(E)dE=C\ exp(-f/T)df 



arctan —\E- Ey2 1 1 - — ^ 



■avctan^lE-Ejl+^^'^ 



(1) 



where F represents the natural width and t the lifetime 
of the depopulated nuclear state. The calculation of v{t) 



was first solved within the so-called Mean Free Path 
Approach (MFPA) [4]. Later also other approaches like 
Molecular Dynamics simulations [5] or Restricted 
Molecular Dynamics [8] have been applied. 

If there are no other sources of additional broadening, 
Eq. (1) and a known function v{t) would allow 
extraction of a nuclear state lifetime and a direct com- 
parison with the experiment. However, there exist two 
further factors which must be taken into account: i) the 
instrumental response function and ii) the Doppler- 
broadening due to thermal motion. Both of them can be 
determined experimentally. In fact, with the double flat 
crystal spectrometers GAMS4 and GAMS5 the determi- 
nation of the instrument response function is possible 
with very high accuracy in the so-called nondispersive 
mode (for more details see Ref. [2]. The thermal 
Doppler-broadening can be determined by scanning a y 
transition, which depopulates a very long lived state 
(t > 10 ps). In this case, it can be assured that the state 
is mostly depopulated, when any recoil motion is slowed 
down and the atom moves with typical thermal 
velocities. Comparing such a scan with a nondispersive 
scan of the same transition (yielding the instrument 
response function) a direct extraction of the thermal 
Doppler-broadening is possible. The knowledge of the 
thermal Doppler broadening allows extraction of a 
target temperature T, which will enter into the calcula- 
tion of the initial and final values of the function v{t) 
within the different slowing down models. Conse- 
quently, the Doppler-broadened lineshape loiE) will 
depend on the target temperature. The total lineshape, 
which needs to be compared to an experiment can be 
written as 



hiE) = UE) * RiE) . 



(2) 
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where hiE) denotes the total Hneshape obtained as a 
convolution of a Doppler-broadened lineshape Id(E) and 
the instrumental response function RiE). Assuming 
the knowledge of v{t) it has been shown that it is possi- 
ble to extract lifetimes of nuclear states in the range of 
10 '^ s to 10 " s using the least squares fitting routine 
GRIDDLE [9]. Also the inverse problem can be consid- 
ered if the lifetime t is known and the function v(t) has 
to be extracted, which allows the study of the atomic 
motion in solids. 

3. Accuracy 

The structure of Eq. (1) suggests that the precision 
with which a nuclear lifetime is extracted depends on 
the accuracy of v(t). Due to the low recoil velocities, the 
theoretical description of the atomic motion can be re- 
stricted to classical mechanics, i.e., to the solution of 
Newton's equations of motion for point mass particles 
interacting via a classical interatomic potential. 

Within the MFPA the slowing down process is 
considered to be on average a sequence of binary atomic 
collisions, where the recoiling atom looses in every 
collision half of its kinetic energy. Between collisions 
the atoms are on average assumed to move with constant 
velocity. The result of this assumption can be summa- 
rized as a step function 



1 



Tk 



14+1(4+1) — /— Vkih), h — / 1 

V2 i=i H-i 



(3) 



where t^ describes the mean free paths between two 
collisions. They can be calculated knowing the density 
of the target material and assuming an atomic interac- 
tion potential [2]. 

A more precise approach for the calculation has been 
introduced to the GRID technique by Kuronen et al. [5], 
who used Molecular Dynamics (MD) simulations for 
the calculation of v{t). This technique is quite common 
in the computer simulation of ion beam solid interaction 
[1 1] and is exact within classical mechanics. The motion 
of the recoiling atom is considered within a simulation 
cell containing A^ atoms, where N=5QQ to 1000. 
Newton's equations of motion of all A' atoms, interacting 
via an interatomic potential, are simultaneously inte- 
grated. However, the computational procedure is rather 
complicated as N{N - l)/2 equations have to be solved. 
Therefore, a slightly simplified approach (Restricted 
Molecular Simulations) has been introduced to the 
GRID-technique [8]. It can be argued that the interac- 
tion of atoms far away from the recoiling nucleus is 
negligible. This allows restriction of the integration to 
only those equations which describe the interaction of 



atoms within a certain sphere surrounding the recoiling 
atom. 

In both MD and RMD simulations, a large number of 
recoil events (typically several thousands) is calculated 
and an average velocity v{t) deduced afterwards. As 
these two models do not directly calculate the average 
velocity they can be used to test the model of an 
isotropic slowing down. This can be done by replacing 
Eq. (1) by the following expression 



UE) = Ci^ d0sin 

7= 1 •'() 



dt 



exp(- t/r) 



[E-E^2{ 1 + c-'vjit)co&e}Y + (rny 



(4) 



Here Ci is a normalisation constant, Vj(t) the velocity 
of the J -th recoiling atom and 6 is the angle under which 
the j -th particle moves with respect to the direction of 
observation. The summation goes over M recoil events. 
The first integration takes account for the isotropic 
alignment of the microcrystals and the second integra- 
tion introduces the emission probability of the second 7 
ray. Comparing lineshapes, which were calculated via 
Eq. (4) to lineshapes calculated via Eq. (1) has shown 
that the isotropic model can be used for polycrystalline 
targets. 

In all three models mentioned above, the basic input 
is defined by the assumption of an interatomic potential. 
Essentially two potentials are used: The MFPA is based 
on a Born-Mayer type repulsive potential [7]. The MD 
and RMD simulations allow implementation of any type 
of potential, but mainly screened Coulomb type poten- 
tials are used [11]. There exists only a sparse knowledge 
on the atomic interaction in solids at kinetic energies of 
several hundreds of eV. In Fig. 2, the function v(t) was 
calculated for ^""Cl atoms recoiling in NaCl, using differ- 
ent theoretical approaches. Within a time interval of up 
to 2 X 10 " s the choice of the interatomic potential 
seems to introduce the dominating uncertainty in the 
slowing down description. Later, the difference of the 
slowing down models becomes dominating. 

In order to quantify the influence of the uncertainty 
in the slowing down description on the lifetime deter- 
mination up to 5 X 10 ''' s, the function v{t) was approx- 
imated by an exponential behavior v{t) = ub exp(-f /T). 
The parameters Vo and T represent the initial recoil 
velocity and the average slowing down time, 
respectively. Although this model is rather simple, it 
gives a surprisingly good fit to the simulated v(t) 
curves for time intervals up to 2 X 10 '^ s. Using this 
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Fig. 2. Time dependence of the average velocity of recoiling Cl-atoms in a NaCl target calculated with three different slowing down models. 
The energies in every plot correspond to the energy of the 7 ray, which initiated the recoil. The thick dashed lines shows the result from Molecular 
Dynamics simulations, the full line the Mean Free Path Approach, and the dash-dotted line the Restricted Molecular Dynamics simulations. In all 
three of them the Born-Mayer potential has been used to model the repulsive atomic interaction. The thick full line describes the results obtained 
from Restricted Molecular Dynamics simulations using a screened Coulomb potential (KrC). 
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approximation the function loiE) can be analytically 
calculated to be 



l{E)dE 



C 



T/t-1 



1 V'-i 


Vo - 



-E-y2 



dE, 



(5) 



where C denotes a normalization constant. The nuclear 
state lifetime enters into equation [5] only as t/T. From 
this it can be concluded that any relative uncertainty of 
the average slowing down time T gives the same relative 
uncertainty for the lifetime extraction. An analysis of 
the calculations presented in Fig. 2 yields an average 
variance of T and therefore t of about 36 %. A similar 
result was obtained in [2], where a detailed comparison 
of several simulations has been performed. According to 
the calculations presented in Fig. 2 it seems that the 
validity of the interatomic potential is the dominating 
uncertainty. However, all comparisons so far have been 
performed for lifetimes up to 10 '^ s and a check for 
longer lifetimes is still missing. This is mainly caused by 
the bad numerical efficiency of MD simulations. 

Up to now considerations have been based on the 
assumption of a fully primarily fed nuclear level of 
interest. In the more general (and more frequent) case, 
the level of interest is fed by several y cascades passing 
several intermediate levels. Here, it is difficult to 
estimate the velocity distribution of the recoiling 
nucleus, when the level of interest is reached. If the 
intensities of all 7 transitions and the lifetimes of all 
intermediate levels are known, this distribution can be 
calculated. In most of the cases this does not apply and 
the uncertainty of this missing knowledge has to be 
quantified with respect to the extraction of nuclear state 
lifetimes from Doppler-broadened lineshapes. For this 
purpose several approaches are used: 

A rather conservative estimation can be done within 
the so-called extreme feeding assumptions. They are 
based on the comparison of the summary intensities 
populating and depopulating the level of interest. If the 
summary intensity of all depopulating transitions is 
greater then the populating one, the missing populating 
intensity can be used to formulate the extreme feeding 
assumptions. In one extreme all missing intensity is 
assigned to a primary feeding and in the other one it 
may be assigned to a transition depopulating a very long 
lived state, which has an energy of about 1 MeV higher 
then the level of interest. In the first case the velocity 
distribution contains a maximum possible number of 
high energetic recoils. This means that in order to obtain 
a certain given Doppler broadening, the lifetime of the 



level of interest has to be rather long. In such a way an 
upper limit for the lifetime estimation can be obtained. 
In the second case, the velocity distribution contains a 
maximum number of low energetic recoils, yielding 
such the lower limit for the lifetime. 

Another possibility for the assignment of the missing 
intensity is given by a so-called x^ analysis of the ex- 
perimental data. Here, a certain parametrization of the 
missing part of the recoil velocity distribution is chosen. 
The parameter values are then included in the fit proce- 
dure for the evaluation of the experimental data. A 
typical parametrization is given by a two step cascade 
(see Fig. 3). Here, all missing populating intensity is 
assigned to a two step-cascade above the level of inter- 
est and the lifetime and the energy of the intermediate 
level are used as parameters. For illustration the two 
step cascade parametrization has been applied to the 
evaluation of data taken with the 1327 keV transition in 
'^'^Gd (Fig. 3). As the evaluation is based on the increase 
of parameters within a least square fit, the data have to 
yield sufficient statistical significance to allow a 
constructive optimization of the recoil velocity distribu- 
tion. However, any result obtained within this approach 
will fall within the limits of the extreme feeding 
assumptions. 

4, Sensitivity 

The simple model of an exponential slowing down 
can be extended to include thermal motion. This can be 
done by folding Eq. (5) additionally with a Doppler- 
broadened lineshape, as it results from thermal motion 
at temperature T with a Maxwell-Boltzmann distribu- 
tion. This takes account for the variance of the initial 
recoil velocity and the final convergence of the slowing 
down to an average thermal velocity. The result ~Id{E, t) 
can not be given in analytical form. However, numeri- 
cally it is easy to handle and can be used to demonstrate 
the sensitivity of the GRID technique to the different 
parameters, which have to be used for an evaluation of 
a measured lineshape. As the evaluation of the data is 
done by means of a least square fitting, the following 
function shall be considered 



^' = 1 



UE,,t+At) UE,t) 



1 



VUE, t) 



^/d(£,, t) 



1 V/d(£,, t) 



(6) 
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Fig. 3. The conception of a two step cascade analysis is demonstrated for a lifetime measurement in ""Gd. The unknown feeding of the 1407 keV 
state is replaced by a two step cascade, where the energy and the lifetime of the intermediate state are varied. For every realisation of this cascade 
a lineshape fit has been performed and the lifetime and the x^ ^e plotted. A most likely feeding scenario could be identified by the minimum x^ 
value. 



Here, £, are discrete points with a certain spacing of AE 
on the energy axis representing the analog to experimen- 
tal data sampling within a scan. The function ^^ has the 
same analytical structure as the partial derivative of the 
optimization function x^ in ^ Is^st square fit. Thus, an 
increase in the values of ^^ means an increase in sensi- 
tivity in lifetime extraction. In Fig. 4 the dependence of 
^^ has been calculated as function of t for assumed 
parameters of Uo = 0.5 A/fs for the initial velocity, 
At/t= 0.1 for a desired precision of the lifetime mea- 
surement, r = 50 fs for the average slowing down time, 
AEi = 10 eV for the scan spacing and Vj = 0.03 A/fs for 
the thermal velocity and auw = 20 eV as FWHM of the 
instrument response function. In general it can be out- 
lined that the highest sensitivity of the technique is given 
in the lifetime region of 10 ''^ to 10 '^ s. At lower life- 
time values the decreasing sensitivity will be compen- 
sated by an increasing influence of the natural width of 
the intermediate level. For lifetimes longer than one 
picosecond the technique loses strongly in sensitivity. 
This can partly be compensated by increased statistics. 
For example, for the extraction of a lifetime of 
T~ 2 X 10 '^ there have to be taken twice as many scans 
as for T» 1 X lO'l 

Using the function ^^ the influence of the instrument 
resolution, the target temperature and the slowing down 



time on the sensitivity have been investigated. Following 
the upper left plot in Fig. 4 the instrument resolution 
shows no substantial influence within the typical ranges 
of variation. The thermal velocity i>r (upper right plot) 
appears to be a strongly limiting factor for the sensitiv- 
ity to lifetimes longer than 10"'^ s, as with increasing 
target temperature more and more statistics is required 
to extract a certain lifetime. The scan interval (lower left 
plot) appears to be of no significance for the sensitivity. 
A strong influence comes from the slowing down time. 
As expected from Eq. (2) an increase in the slowing 
down time would allow to extend the technique towards 
longer lifetimes. 



5. Conclusions 

The GRID technique has been shown to be a sensi- 
tive experimental technique for the investigation of nu- 
clear state lifetimes in the range of 10 '^ to 10 '^ s. It 
requires the knowledge of atomic motion with several 
hundreds of eV kinetic energy and the feeding scenario 
of the nuclear level of interest. 

The relative uncertainty of the theoretical description 
of atomic motion can be quantified in certain cases. It 
results mainly from the unknown atomic interaction. 
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Fig. 4. The function ^^ is shown for different parameter sets for the calculation oilDiE;, t). In every plot one parameter has been varied and 
compared to a reference calculation (dashed line), which was calculated for £^=1 MeV, n; = 0.5 X 10^ m/s, At/t=0.1, 7=50 X 10"'^ s, 
AE, = 10 eV, (7hw = 20 eV and vp = 0.03 X 10' m/s. 



Lifetimes in the range of 10 '^ to 10 " s measured with 
the GRID technique have the same relative uncertainty 
as the slowing down description. Recently, efforts have 
been made to improve the current status and to study 
explicitly the atomic interaction using the GRID tech- 
nique itself. This can be done using single crystal targets 
(Crystal-GRID), which allows a significant increase in 
the amount of information obtained. Further details can 
be found in Refs. [12, 13]. 

An additional uncertainty for the lifetime extraction is 
caused by an unknown feeding scenario for the level of 
interest. It can be quantified using the concept of 
extreme feeding scenarios, which yields the maximum 
variation of the measured lifetime value as it can be 
produced by the unknown feeding. Using the approxi- 
mation of a two step cascade for the missing feeding one 
can find a most likely recoil velocity distribution and 
extract the most likely directly from the experimental 
data. Alternativelly, the missing feeding can be calcu- 
lated using statistical models of the nucleus [14]. 
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